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Abstract. Early 17th-century mathematical publications of Johann Faulhaber 
contain some remarkable theorems, such as the fact that the r-fold summation 
of l m , 2 m , . . . , n m is a polynomial in n{n + r) when m is a positive odd number. 
The present paper explores a computation-based approach by which Faulhaber 
may well have discovered such results, and solves a 360-year-old riddle that Faul- 
haber presented to his readers. It also shows that similar results hold when we 
express the sums in terms of central factorial powers instead of ordinary powers. 
Faulhaber's coefficients can moreover be generalized to factorial powers of non- 
integer exponents, obtaining asymptotic series for 1° + 2 a + ■ ■ • + n a in powers 
of n _1 (n + 
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Johann Faulhaber of Ulm (1580-1635), founder of a school for engineers early in the 17th 
century, loved numbers. His passion for arithmetic and algebra led him to devote a consid- 
erable portion of his life to the computation of formulas for the sums of powers, significantly 
extending all previously known results. He may well have carried out more computing than 
anybody else in Europe during the first half of the 17th century. His greatest mathematical 
achievements appear in a booklet entitled Academia Algebrae (written in German in spite 
of its latin title), published in Augsburg, 1631 [2]. Here we find, for example, the following 
formulas for sums of odd powers: 

l 1 + 2 1 H + n 1 = N , N = (n 2 + n)/2 ; 

l 3 + 2 3 + • • • + n 3 = N 2 ; 

l 5 +2 5 + --- + n 5 = (4N 3 -N 2 )/3; 

l 7 + 2 7 + • • • + n 7 = (12iV 4 - 8N 3 + 2N 2 ) /6 ; 

l 9 + 2 9 + • • • + n 9 = (16N 5 - 2(W 4 + 12N 3 - 3iV 2 )/5 ; 
l 11 + 2 11 + • • • + n 11 = (32iV 6 - 64iV 5 + 68N 4 - AON 3 + 5iV 2 ) /6 ; 
l 13 + 2 13 H + n 13 = (960N 7 - 2800iV 6 + 4592N 5 - 4720iV 4 + 27647V 3 

- 691N 2 )/105 ; 

l 15 + 2 15 H h n 15 = (192N 8 - 768iV 7 + 1792N 6 - 2816iV 5 + 2872A^ 4 

- 1680A^ 3 + 42(W 2 )/12; 

l 17 + 2 17 H + n 17 = (1280A^ 9 - 6720iV 8 + 21120A^ 7 - 46880iV 6 + 72912A^ 5 

- 74220A^ 4 + 43404iV 3 - 10851A^ 2 ) /45 . 

Other mathematicians had studied En 1 , En 2 , . . . , En 7 and he had previously gotten as far 
as En 12 ; but the sums had always previously been expressed as polynomials in n, not N. 

Faulhaber begins his book by simply stating these novel formulas and proceeding to 
expand them into the corresponding polynomials in n. Then he verifies the results when 
n = 4, iV = 10. But he gives no clues about how he derived the expressions; he states only 
that the leading coefficient in En 2m_1 will be 2 m_1 /m, and that the trailing coefficients 
will have the form 4a m iV 3 — a m N 2 when m > 3. 

Faulhaber believed that similar polynomials in N, with alternating signs, would con- 
tinue to exist for all m, but he may not really have known how to prove such a theorem. In 
his day, mathematics was treated like all other sciences; it was sufficient to present a large 
body of evidence for an observed phenomenon. A rigorous proof of Faulhaber's assertion 
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was first published by Jacobi in 1834 [6]. A. W. F. Edwards showed recently how to obtain 
the coefficients by matrix inversion [1], based on another proof given by L. Tits in 1923 [8]. 
But none of these proofs use methods that are very close to those known in 1631. 

Faulhaber went on to consider sums of sums. Let us write E r n m for the r-fold sum- 
mation of mth powers from 1 to n; thus, 

E°n m = n m ; E r+1 n m = E r l m + E r 2 m + • • • + E r n m . 

He discovered that E r n 2m can be written as a polynomial in the quantity 

N r = (n 2 + rn)/2, 

times E r n 2 . For example, he gave the formulas 



E 2 n 4 


= (4N 2 - 1) E 2 n 2 /5 ; 


E 3 n 4 


= (4iV 3 - 1) E 3 n 2 /7 ; 


E 4 n 4 


= (6iV 4 - 1) E 4 ?i 2 /14 ; 


E 6 n 4 


= (4iV 6 + 1) E 6 n 2 /15 ; 


E 2 n 6 


= (6iV 2 2 - 5N 2 + 1) E 2 n 2 /7 ; 


E 3 n 6 


= (lOiVf - 10N 3 + 1) E 3 n 2 /21 ; 


E 4 n 6 


= (4iVf - 4AT 4 - 1) E 4 n 2 /14 ; 


E 2 n 8 


= (16AT 3 - 28AT 2 + 18JV 2 - 3) E 2 n 2 /15 



He also gave similar formulas for odd exponents, factoring out E r n 1 instead of E r n 2 : 

E 2 n 5 = (8AT| - 2N 2 - 1)EV/14 ; 

E 2 n 7 = (40iV 2 3 - 40AT| + 6N 2 + 6)E 2 nV60 . 

And he claimed that, in general, E r n m can be expressed as a polynomial in A^ r times either 
E r n 2 or E 7 "?! 1 , depending on whether m is even or odd. 

Faulhaber had probably verified this remarkable theorem in many cases including 
E n n 6 , because he exhibited a polynomial in n for E n n 6 that would have been quite 
difficult to obtain by repeated summation. His polynomial, which has the form 

6n 17 + 561n 16 + • • • + 1021675563656n 5 + 96598656000n 

2964061900800 ' 

turns out to be absolutely correct, according to calculations with a modern computer. 
(The denominator is 171/120. One cannot help thinking that nobody has ever checked 
these numbers since Faulhaber himself wrote them down, until today.) 
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Did he, however, know how to prove his claim, in the sense that 20th century math- 
ematicians would regard his argument as conclusive? He may in fact have known how to 
do so, because there is an extremely simple way to verify the result using only methods 
that he would have found natural. 

Reflective functions. Let us begin by studying an elementary property of integer func- 
tions. We will say that the function f(x) is r-reflective if 

fix) = f(y) whenever x + y + r = ; 

and it is anti-r-reflective if 

f(x) = —f(y) whenever x + y + r = 0. 

The values of x, y, r will be assumed to be integers for simplicity. When r = 0, reflective 
functions are even, and anti-reflective functions are odd. Notice that r-reflective functions 
are closed under addition and multiplication; the product of two anti-r-reflective functions 
is r-reflective. 

Given a function /, we define its backward difference V/ in the usual way: 

V/(») = f(x) - f{x - 1) . 
It is now easy to verify a simple basic fact. 

Lemma 1. If / is r-reflective then V/ is anti-(r — 1) -reflective. If f is anti-r-reflective 
then V/ is (r — 1) -reflective. 

Proof. If x + y + (r — 1) =0 then x + iy — 1) + r = and (x — 1) + y + r = 0. Thus 
fix) — ±/(y — 1) and fix — 1) = ±/(y) when / is r-reflective or anti-r-reflective. □ 

Faulhaber almost certainly knew this lemma, because [2, folio D.iii recto] presents a 
table of n 8 , Vn 8 , . . . , V 8 n 8 in which the reflection phenomenon is clearly apparent. He 
states that he has constructed "grosse Tafeln," but that this example should be "alles 
gnugsam vor Augen sehen und auf hohere quantiteten [exponents] continuiren konde." 

The converse of Lemma 1 is also true, if we are careful. Let us define E as an inverse 
to the V operator: 

nU) \C-f(0) /(n+1), ifn<0. 

Here C is an unspecified constant, which we will choose later; whatever its value, we have 

VE/(n) = S/(n) - E/(n - 1) = f(n) 

for all n. 
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Lemma 2. If f is r-reflective, there is a unique C such that £/ is anti-(r + 1) -reflective. 
If f is anti-r-reflective, then £/ is (r + 1) -reflective for all C. 

Proof. If r is odd, £/ can be anti-(r + l)-refiective only if C is chosen so that we have 
^f{~~ ( r + =0. If r is even, £/ can be anti-(r + l)-reflective only if E/(— r/2) = 

-£/(-r/2 - 1) = _(E/(-r/2) - /(-r/2)); i.e., £/(-r/2) = J/(-r/2). 

Once we have found a; and y such that x + y + r + l = and E/(x) = — E/(y), it 
is easy to see that we will also have E/(x — 1) = —Hf(y + 1), if / is r-reflective, since 
£/0r) - E/(x - 1) = /(x) = /(» + 1) = E/fo + 1) - E/(y). 

Suppose on the other hand that / is anti-r- reflective. If r is odd, clearly E/(x) = E/(y) 
if x = y = — (r + l)/2. If r is even, then /(—r/2) = 0; so E/(x) = E/(y) when a; = —r/2 
and y = —r/2 — 1. Once we have found x and y such that x + y + r + 1 = and 
E/(x) = E/(y), it is easy to verify as above that E/(x — 1) = E/(y + 1). □ 

Lemma 3. If f is any even function with /(0) = 0, the r-fold repeated sum E r / is r- 
reflective for all even r and anti-r -reflective for all odd r, if we choose the constant C = 
in each summation. If f is any odd function, the r-fold repeated sum E r / is r-reflective 
for all odd r and anti-r -reflective for all even r, if we choose the constant C = in each 
summation. 

Proof. Note that /(0) = if / is odd. If /(0) =0 and if we always choose C = 0, it is 
easy to verify by induction on r that T, r f(x) = for — r < x < 0. Therefore the choice 
C = always agrees with the unique choice stipulated in the proof of Lemma 2, whenever 
a specific value of C is necessary in that lemma. □ 

When m is a positive integer, the function f(x) = x m obviously satisfies the condition 
of Lemma 3. Therefore we have proved that each function E r n m is either r-reflective or 
anti-r-reflective, for all r > and m > 0. And Faulhaber presumably knew this too. 
His theorem can now be proved if we supply one small additional fact, specializing from 
arbitrary functions to polynomials: 

Lemma 4. A polynomial f(x) is r-reflective if and only if it can be written as a polynomial 
in x(x + r); it is anti-r-reflective if and only if it can be written as (x + r/2) times a 
polynomial in x(x + r). 

Proof. The second statement follows from the first, because we have already observed 
that an anti-r-reflective function must have /(—r/2) = and because the function x + r/2 
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is obviously anti-r-reflective. Furthermore, any polynomial in x(x + r) is r-reflective, 
because x(x + r) = y(y + r) when x + y + r = 0. Conversely, if f(x) is r-reflective we 
have f(x — r/2) = f(—x — r/2), so g(x) = f(x — r/2) is an even function of x; hence 
g(x) = h(x 2 ) for some polynomial h. Then f(x) = g(x + r/2) = h(x(x + r) + r 2 /4) is a 
polynomial inx(x + r). □ 

Theorem (Faulhaber). There exist polynomials g r ^ m for all positive integers r and m 
such that 

E r n 2m-i = gr 2m+i ( n{n + r) j S r n i ^ s V m = </ r>2m (n(n + r)) E r n 2 . 

Proof. Lemma 3 tells us that E r n m is r-reflective if m + r is even and anti-r-reflective if 
m + r is odd. 

Note that E 7 '?! 1 = ("^[) • Therefore a polynomial in n is a multiple of E r n 1 if and 
only if it vanishes at — r, ...,—1,0. We have shown in the proof of Lemma 3 that E r n m 
has this property for all m; therefore E r n m /E r n 1 is an r-reflective polynomial when m is 
odd, an anti-r-reflective polynomial when m is even. In the former case, we are done, by 
Lemma 4. In the latter case, Lemma 4 establishes the existence of a polynomial g such 
that E r n m /E r n 1 = (n + r/2)g(n(n + r)). Again, we are done, because the identity 

r + 2 

is readily verified. □ 

A plausible derivation. Faulhaber probably didn't think about r-reflective and anti-r- 
reflective functions in exactly the way we have described them, but his book [2] certainly 
indicates that he was quite familiar with the territory encompassed by that theory. 

In fact, he could have found his formulas for power sums without knowing the theory 
in detail. A simple approach, illustrated here for En 13 , would suffice: Suppose 

14E?i 13 = n 7 (n + l) 7 - S(n) , 

where S(n) is a 1-reflective function to be determined. Then 

14n 13 = n 7 (n + l) 7 - (n - 1) V - VS(n) 

= 14n 13 + 7072 11 + 42n 9 + 2n 7 - VS(n) , 
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and we have 

S(n) = 70EU 11 + 42n 9 + 2En 7 . 

In other words 

En 13 = ^ N 7 _ 5En ll _ 3En 9 _ I En 7 

7 7 
and we can complete the calculation by subtracting multiples of previously computed 
results. 

The great advantage of using polynomials in N rather than n is that the new formulas 
are considerably shorter. The method Faulhaber and others had used before making this 
discovery was most likely equivalent to the laborious calculation 

En 13 = _L_ n 14 + ±3^12 _ 26En ll + l|3 En 10 _ 143En 9 + 429^8 

+ i^En 7 + ifEn 6 - 143En 5 + if En 4 - 26En 3 + f En 2 - En 1 + ±n ; 

the coefficients here are ^(J^K -n(ii)> • • • ' TzCo)- 

To handle sums of even exponents, Faulhaber knew that 

En 2m = ^_J2^ + ^2 + . . . + amN m) 

2m + 1 

holds if and only if 

En 2m+1 = Ol n2 + 02 N s + + _On^ N m+l 

2 3 m + 1 

Therefore he could get two sums for the price of one [2, folios Civ verso and D.i recto]. It is 
not difficult to prove this relation by establishing an isomorphism between the calculations 
of En 2m+1 and the calculations of the quantities S 2m = ((2m + l)En 2m )/ (n + |); for 
example, the recurrence for En 13 above corresponds to the formula 

S 12 = 64iV 6 - 5S 10 - 3S 8 - ±S 6 , 

which can be derived in essentially the same way. Since the recurrences are essentially 
identical, we obtain a correct formula for En 2m+1 from the formula for S2m if we replace N k 
by N k+1 /(k+ 1). 

Faulhaber's cryptomath. Mathematicians of Faulhaber's day tended to conceal their 
methods and hide results in secret code. Faulhaber ends his book [2] with a curious exercise 
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of this kind, evidently intended to prove to posterity that he had in fact computed the 
formulas for sums of powers as far as En 25 although he published the results only up to 
En 17 . 

His puzzle can be translated into modern notation as follows: Let 



9 8 ai 7 n H h a 2 n z + a\n 



d 

where the a's are integers having no common factor and d = au + • ■ • + a 2 + a\. Let 

25 _ A 26 n 26 + ■■■ + A 2 n 2 + Am 
En - - 

be the analogous formula for En 25 . Let 

v 22 (b 10 N 10 -b 9 N 9 + --- + b ) 2 

En = - - - En , 

bio -o 9 -\ \-t>o 

En 23 = (Clonl °- C9iV9 + --- + Co) En 3 , 
cio - c 9 H h c 

a," = (*>""-"^"' + --'fo) a ,. , 

dn - dio H h do 

En- = (^n"-e 10 iV» + ...-e ) Eti3 
en — eio H h e 

where the integers b k , c k , d k , e k are as small as possible so that b k ,c k ,d k ,e k are mul- 
tiples of 2 k . (He wants them to be multiples of 2 k so that b k N k , c k N k , d k N k , e k N h 
are polynomials in n with integer coefficients; that is why he wrote, for example, En 7 = 
(12N 2 -8N + 2)N 2 /6 instead of (6N 2 - AN + l)N 2 /3. See [2, folio D.i verso].) Then 
compute 

Xl = (c 3 - ai 2 )/7924252; 
X2 = (6 5 + aio) /112499648 ; 
x 3 = (an -bg- ci)/2945002 ; 
x 4 = (a M + c 7 )/120964; 

x 5 = (A 26 a 11 - D + a 13 + d u + e n )/199444 . 

These values (xi, x 2 , £3, X4, x 5 ) specify the five letters of a "hochgeriihmte Nam," if we use 
five designated alphabets [2, folio F.i recto]. 
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It is doubtful whether anybody solved this puzzle during the first 360 years after its 

publication, but the task is relatively easy with modern computers. We have 

aio = 532797408, an = 104421616, a 12 = 14869764, a 13 = 1526532, a M = 110160; 

6 5 = 29700832 , b 9 = 140800 ; 

ci = 205083120 , c 3 = 344752128 , c 7 = 9236480 ; 

dn = 559104; en = 86016; A 26 = 42 ; D = 1092 . 

The fact that x 2 = (29700832 + 532797408)/112499648 = 5 is an integer is reassuring: We 
must be on the right track! But alas, the other values are not integral. 

A bit of experimentation soon reveals that we do obtain good results if we divide all 
the c k by 4. Then, for example, x x = (344752128/4 - 14869764)/7924252 = 9, and we also 
find X3 = 18, £4 = 20. It appears that Faulhaber calculated E 9 n 8 and En 22 correctly, and 
that he also had a correct expression for En 23 as a polynomial in N; but he probably never 
went on to express En 23 as a polynomial in n, because he would then have multiplied his 
coefficients by 4 in order to compute c 6 iV 6 with integer coefficients. 

The values of (xi, x 2 , £3, X4) correspond to the letters IESU, so the concealed name 
in Faulhaber's riddle is undoubtedly IESUS (Jesus). 

But his formula for x$ does not check out at all; it is way out of range and not an 
integer. This is the only formula that relates to En 24 and En 25 , and it involves only the 
simplest elements of those sums — the leading coefficients A 2 q, D, dn, en. Therefore we 
have no evidence that Faulhaber's calculations beyond En 23 were reliable. It is tempting 
to imagine that he meant to say 'A 26 cin/D' instead of 'A 26 an — D' in his formula for x$, 
but even then major corrections are needed to the other terms and it is unclear what he 
intended. 

All-integer formulas. Faulhaber's theorem allows us to express the power sum En m 
in terms of about |m coefficients. The elementary theory above also suggests another 
approach that produces a similar effect: We can write, for example, 

»=(?); 

« 3 = 6(?) + C); 

n 5 = 120 ( n J 2 ) +30( n J 1 ) + (?) ; 
(It is easy to see that any odd function g(n) of the integer n can be expressed uniquely as 
a linear combination 

g(n) = ai Q) + as^t 1 ) + as( n+ 5 2 ) + ■ ■ ■ 
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of the odd functions (™), ( n + 1 ), ( n J 2 ), because we can determine the coefficients 
ai, 03, 05, . . . successively by plugging in the values n = 1, 2, 3, . . . . The coefficients will 
be integers iff g(n) is an integer for all n.) Once g(n) has been expressed in this way, we 
clearly have 

Zg(n) = 01 ( n J x ) + a 3 ( n t 2 ) + a 5 ( n + 3 ) + • • • . 
This approach therefore yields the following identities for sums of odd powers: 

En* =rr); 

En 3 =6(»+ 2 ) + (»+ 1 ); 

En 5 = 120 ( n + 3 ) + 30 ( n + 2 ) + ( n + 1 ) ; 

En 7 = 5040 ( n J 4 ) + 1680 ( n J 3 ) + 126 ( n + 2 ) + ("J 1 ) ; 

En 9 = 362880 ( n + 5 ) + 151200 ("J 4 ) + 17640 ("J 3 ) + 510 ( n + 2 ) + ("J 1 ) ; 

En 11 = 39916800 ( n + 6 ) + 19958400 ( n + 5 ) + 3160080 ( n J 4 ) + 168960 ( n J 3 ) 

+ 2046( n + 2 ) + ( n + 1 ); 
En 13 = 6227020800 ( n + 7 ) + 3632428800 ( n + 6 ) + 726485760 ( n + 5 ) 

+ 57657600 ( n + 4 ) + 1561560 ( n + 3 ) + 8190 ( n + 2 ) + ( n + 1 ) . 
And repeated sums are equally easy; we have 

E V = (»£) , E^ 3 = 6 + QZ) > etc. 

The coefficients in these formulas are related to what Riordan [R, page 213] has called 
central factorial numbers of the second kind. In his notation 



E T ( m > k ) x[k] » xlk] = x ( x + i ~ 1 )( x + i ~ 2 ) • • • ( x + i + !) 

k=i 

when m > 0, and T(m, fc) = when m — is odd; hence 



n 2 - 1 = £ (2* - 1)! T(2m, 2fe) + * 

fc=i ^ ' 

En 2 ™" 1 = ^ (2Jfe - 1)! T(2m, 2fc) T 
fc=i ^ 



n + /c 
2Jfe 

The coefficients T(2m, 2/c) are always integers, because the identity x^ h+2 ^ = x^ k \x 2 — k 2 /4) 
implies the recurrence 

T(2m + 2, 2k) = k 2 T(2m, 2k) + T(2m, 2k -2). 
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The generating function for these numbers turns out to be 



oo / m 



cosh(2xsinh(y/2)) = ^ I ^T(2m,2/c)x 



2k\ y 



2m 



(2m)! 



m=0 y k=0 

Notice that the power-sum formulas obtained in this way are more "efficient" than 
the well-known formulas based on Stirling numbers (see [5, (6.12)]): 



En r 



ml (n + 1 
k J U + l 



5> 



k fn + k 
k+1 



The latter formulas give, for example, 

En 7 = 5040 ( n + x ) + 15120 ^J 1 ) + 16800 ( n + x ) + 8400 ("J 1 ) + 1806 

+ 126P+ 1 ) + 

= 5040 ( n J 7 ) - 15120 ( n + 6 ) + 16800 ( n + 5 ) - 8400 ( n J 4 ) + 1806 ( n J 3 ) 
- 126(»f) + C^) • 

There are about twice as many terms, and the coefficients are larger. (The Faulhaberian 
expression En 7 = (6iV 4 - 4iV 3 + N 2 )/3 is, of course, better yet.) 

Similar formulas for even powers can be obtained as follows. We have 



n 2 = nQ) 

n 4 _ (n+l 



n 



6n( n + 1 ) +n( n 1 ) = 12U 2 (n) + U^n) , 

120n( n J 2 ) + SOn^ 1 ) +nQ)= 360*7 3 (ra) + 60U 2 (n) + U x {n) , 



etc., where 



U k (n) = 



n (n + k — 1 



k V 2k -1 



n + k 
2k 



+ 



n + k — 1 
2k 



Hence 



En 2 
En 4 
En 6 
En 8 



Ti(n), 

12T 2 (ra) + Ti(n) , 
360T 3 (n) + 60T 2 + T 1 (n), 
20160T 4 (n) + 5040T 3 (n) + 252T 2 (n) + Ti(n) , 
En 10 = 1814400T 5 (n) + 604800T 4 (n) + 52920T 3 (n) + 1020T 2 (n) + T^n) , 
En 12 = 239500800T 6 (n) +99792000T 5 (n) + 12640320T 4 (n) 
+ 506880T 3 (n) + 4092T 2 (n) + Ti(ra) , 
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etc., where 

(n + k + l\ /n + k\_2n + l/n + k\ 
V 2k +1 ) + \2k+l) ~ 2k + 1 V 2k )' 
Curiously, we have found a relation here between Sn 2m and En 2m_1 , somewhat anal- 
ogous to Faulhaber's relation between En 2m and En 2m+1 : The formula 

2^1 = H 2 J + H 4 J + --- + ^^ 2m ) 
holds if and only if 

2m-i 3 A^ + A 5 /n + 2\ 2m + 1 /n + m\ 

4. Reflective decomposition. The forms of the expressions in the previous section lead 
naturally to useful representations of arbitrary functions f(n) defined on the integers. It 
is easy to see that any f(n) can be written uniquely in the form 

ft \ sr ( n +l k / 2 \\ 

k>o v fc / 
for some coefficients a^; indeed, we have 

a k =V k f([k/2\). 

(Thus a = /(0), ai = /(0) - /(-l), a 2 = /(l) - 2/(0) + /(-l), etc.) The a fc are integers 
iff f{n) is always an integer. The a k are eventually zero iff / is a polynomial. The d2k are 
all zero iff / is odd. The a,2k+i are all zero iff / is 1-reflective. 
Similarly, there is a unique expansion 

f(n) = b T Q (n) + hUtin) + b 2 T 1 (n) + b 3 U 2 (n) + b 4 T 2 (n) + • • • , 

in which the bk are integers iff f(n) is always an integer. The b 2 k are all zero iff / is even 
and /(0) = 0. The b 2 k+i are all zero iff / is anti-l-reflective. Using the recurrence relations 

VT fc (n) = U k (n), VU k {n) = T fc _i(n - 1) , 

we find 

a k = V k f([k/2\) = 2b k . 1 + (-l) k b k 

and therefore 

k 

6 fc = ^(-l)ri/ 2 l + L*/2J 2 *-i ai . 
i=o 

In particular, when /(n) = 1 for all n, we have b k = (— 1) l k / 2 l2 k . The infinite series is 
finite for each n. 



T k (n) 
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Theorem. If f is any function defined on the integers and if r, s are arbitrary integers, 
we can always express f in the form 

f(n) = g(n) + h(n) 

where g(n) is r-re&ective and h(n) is anti-s-reflective. This representation is unique, except 
when r is even and s is odd; in the latter case the representation is unique if we specify 
the value of g or h at any point. 

Proof. It suffices to consider < r, s < 1, because f(x) is (anti)-r-reflective iff f(x + a) is 
(anti)-(r + 2a)-reflective. 

When r = s = 0, the result is just the well known decomposition of a function into 
even and odd parts, 

g(n) = i(/(n) + f(-n)) , h(n) = \{f{n) - f(-n)) . 

When r = s = 1, we have similarly 

g(n) = \ (f(n) + /(-l - n)) , h(n) = \ (f(n) - f(-l - n)) . 

When r = 1 and s = 0, it is easy to deduce that h(0) = 0, g(0) = /(0), h(l) = 
/(0) - /(-l), g(l) = /(l) - /(0) + /(-l), h{2) = /(l) - /(0) + /(-l) - /(-2), g(2) = 
/(2)-/(l) + /(0)-/(-l) + /(-2),etc. 

And when r = and s = 1, the general solution is g(0) = f(0) — C, h(0) = C, 
g(l) = /(-l) + C, h(l) = /(l) - /(-l) - C, g{2) = /(l) - /(-l) + /(-2) - C, fc(2) = 
/(2)-/(l) + /(-l)-/(-2) + C,etc. □ 

When f(n) = J2k>o a k( n+ ^k^ 2 ^)i the case r =1 and s = corresponds to the decom- 
position 

00 ( n _|_ fc\ 00 
f( n ) = 5Z tt2fc ( 2A; J' M") = $^a2fc+i 

fc=0 ^ ' fc=0 

Similarly, the representation f(n) = J2k>o^ 2k ^ k ( n ) + Sfc>o ^ 2fc + 1 ^ fc + 1 ( n ) corresponds to 
the case r = 0, s = 1, C = f(0). 

Back to Faulhaber's form. Let us now return to representations of En m as polynomials 

in n(n + 1). Setting u = 2N = n 2 + n, we have 
En = \u = lA^u 

Zn s = ±u 2 =\{Af ) u 2 + Af ) u) 

En 5 = \ {u 3 - \u 2 ) = | (4V + 4 3) « 2 + 4 3) «) 

En 7 = i( W 4 - |w 3 + fw 2 ) = |(4 4) w 4 + A< 4 V + Afu 2 + A^u) 
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n + k 
2k + 1 



and so on, for certain coefficients A y k ' . 

Faulhaber never discovered the Bernoulli numbers; i.e., he never realized that a single 
sequence of constants B , Bi, B 2 , ■ ■ ■ would provide a uniform formula 

Sn m = ^r(S n m+1 - ( m + 1 )B 1 n m + ( m + 1 )B 2 n m ~ 1 -■■■ + {-l) m { m ^)B m n) 

for all sums of powers. He never mentioned, for example, the fact that almost half of 
the coefficients turned out to be zero after he had converted his formulas for En m from 
polynomials in N to polynomials in n. (He did notice that the coefficient of n was zero 
when m > 1 was odd.) 

However, we know now that Bernoulli numbers exist, and we know that B3 = B$ = 
B 7 = ■ ■ ■ = 0. This is a strong condition. Indeed, it completely defines the constants 
in the Faulhaber polynomials above, given that 4 m) = 1- 

For example, let's consider the case m = 4, i.e., the formula for En 7 : We need to find 
coefficients a = b = A^\ c = A^ such that the polynomial 

n 4 (n + l) 4 + an 3 (n + l) 3 + bn 2 (n + l) 2 + cn(n + 1) 

has vanishing coefficients of n 5 , n 3 , and n. The polynomial is 

n 8 + An 7 + 6n 6 + 4n 5 + n 4 

+ an 6 + 3an 5 + 3an 4 + an 3 

+ bn 4 + 2bn 3 + n 2 

+ en 2 + cn ; 

so we must have 3a + 4 = 2b + a = c = 0. In general the coefficient of, say, n 2rn ~ 5 in the 
polynomial for 2m£n 2m_1 is easily seen to be 

(-)4 m) + ( m 3 - 1 )A^ m) + ( m " 2 )4 m) • 
Thus the Faulhaber coefficients can be defined by the rules 

4 w) = i; y( 7 W ~ 3 W ) = o, k>o. (*) 

p< Q \2k + 1 - 2jJ J ' v 1 

(The upper parameter will often be called w instead of m, in the sequel, because we will 
want to generalize to noninteger values.) Notice that (*) defines the coefficients for each 



14 



exponent without reference to other exponents; for every integer k > 0, the quantity A { k 
is a certain rational function of w. For example, we have 

-A[ w) =w(w- 2) /6, 

4"° =«;(«; -l)(w-3)(7i<;- 8)/360, 
-4™ } = w(w -l)(w- 2)(w - 4)(31w 2 - 89w + 48)/15120 , 

A[ w) = w(w -l)(w- 2){w -S)(w- 5)(127w 3 - 691w 2 + 1038w - 384)/6048000 , 

and in general A™ is w— = w(w — 1) ... (w — k + 1) times a polynomial of degree /c, with 
leading coefficient equal to (2 — 2 2fc )£?2fc/(2/c)! ; if k > 0, that polynomial vanishes when 
u; = + 1. 

Jacobi mentioned these coefficients 4"^ in [6], although he did not consider the 
recurrence (*), and he tabulated them for m < 6. He observed that the derivative of En m 
with respect to n is mSn m_1 + B m ; this follows because power sums can be expressed in 
terms of Bernoulli polynomials, 

£n™ = ^(B m+1 (n + 1) - B m+i (0)) , 

and because B' m (x) = mB m _ 1 (x). Thus Jacobi obtained a new proof of Faulhaber's 
formulas for even exponents 

En 2 = i(i4 2) n+i4 2 ))(2n + l), 

En 4 = i(|4 3 V + 14 3) ^ + |4 3) )(2n + 1) , 

En 6 = \ (§4V + |4 4 V + |4 4) « + |4 4) ) (2n + 1) , 
etc. (The constant terms are zero, but they are shown explicitly here so that the pattern 
is plain.) Differentiating again gives, e.g., 

En 5 = ( (4 . 3 4V + 3-2 A^u + 2-1 4 4) ) (2n + l) 2 

6 ■ 7 ■ 8 

+ 2(44 4) « 3 + 34V + 24 4) W + l4 4) )) - ±B 6 

= ( 8 ■ 7 4V + (6 • 5 4 4) + 4-3 4 4 V 2 + (4 • 3 4 4) + 3 • 2 ^Ju 

6-7-8 

+ (2.l4 4 >+2.l4 4 >))-I J B 6 . 
This yields Jacobi' s recurrence 

(2w - 2k)(2w -2k- 1)4" ) + (w-k + l)(w - A;)4-i = M 2 ™ - 1)4" _1) , (**) 

15 



which is valid for all integers w > k + 1 so it must be valid for all w. Our derivation of 
(**) also allows us to conclude that 



by considering the constant term of the second derivative of En 2m_1 . 

Recurrence (*) does not define A^\ except as the limit of when w — > m. But 
we can compute this value by setting w = m + 1 and k = m in (**), which reduces to 

2A ( r +i) = (2m + 2)(2m + 1)A M 

because A^ 1 ^ = 0. Thus 

= B 2 m , integer m > . 



Solution to the recurrence. An explicit formula for A^ can be found as follows: We 
have 

En 2 ™" 1 = ^ (B 2m (n + 1) - B 2m ) = ^{A^u™ + ■■■ + A^\u) , 
and n + 1 = (y/l + 4u + l)/2; hence, using the known values of A^\ we obtain 

a closed form in terms of Bernoulli polynomials. (We have used the fact that A^^_ x = 
= ■ • • = 0, together with the identity B n (x + 1) = (— l) n B n (— x) .) Expanding the 
right side in powers of u gives 

l 



using equation (5.70) of [5]. Setting j + I = m — k finally yields 



\m — k — j J \ j J m — k + j 
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This formula, which was first obtained by Gessel and Viennot [4] , makes it easy to confirm 
that A^™\ = and A^} 2 = ( 2 ™)£?2m-2 ? and to derive additional values such as 



A 



(m) 
m— 3 



(?) 



A 



(m) 
m— 4 



2 

2m 



B 2m _ 2 = -2A { ™1 2 , m>3; 



I ,-B 2m _4 + 5l \B 2 m 2 



2m 



m > 4 . 



The author's interest in Faulhaber polynomials was inspired by the work of Ed- 
wards [1], who resurrected Faulhaber's work after it had been long forgotten and underval- 
ued by historians of mathematics. Ira Gessel responded to the same stimulus by submitting 
problem E3204 to the Math Monthly [3] regarding a bivariate generating function for Faul- 
haber's coefficients. Such a function is obtainable from the univariate generating function 
above, using the standard generating function for Bernoulli polynomials: Since 



X + l 



v 2m 



1 



(2m)! 



x + l 



5J + 5 5> 



x + l 



2^ "'V 2 ) m\ 2^ V 2 

2e (x+i)z/2 ^ ze -(*+i)*/2 ^ ZC osh(xz/2) 
2(e* - 1) ~ 2(e"* - 1) ~ 2 sinh(^/2) 



m! 



we have 



k,n 



{m) u m-k z 



2m 



E- 4 

k,w 



(m) k z 
k u 



(2m)! 

2m 



= E* 



2m 



2m 



(2m) 



2 cosh (a/1 + 4«2;/2) 
2 sinh (z/2) 



(2m) 



z^pu cosh (yV+Iz/2) 
2sinh(z v /u/2) 



The numbers A^ are obtainable by inverting a lower triangular matrix, as Edwards 
showed; indeed, recurrence (*) defines such a matrix. Gessel and Viennot [4] observed that 
we can therefore express them in terms of a k x k determinant, 



iH 



(1 — w) ... (k — w) 



^w — k+1^ — 







'w-k+2^ ^w-k+2^ ^w-k+2^ 

(w—l\ (w-l\ ( w ~ 1 '\ 

\2k-l) \2k-2>) \2k-5) 

I w \ I w \ I w \ 

\2k+l) \2k-l) \2k-3) 








(V) 

(3) 
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When w and k are positive integers, Gessel and Viennot proved that this determinant is 
the number of sequences of positive integers aia2a 3 . . . a^k such that 

a 3j _ 2 < a 3j _i < a 3j - < w - k + j , for 1 < j < k , 
a 3 j_ 2 < a 3j+ i , a 3j _i < a 3j+3 , for 1 < j < k . 

In other words, it is the number of ways to put positive integers into a /c-rowed triple 
staircase such as 



with all rows and all columns strictly increasing from left to right and from top to bottom, 
and with all entries in row j at most w — k + j . This provides a surprising combinatorial 
interpretation of the Bernoulli number B 2 m when w = m + 1 and k = m — 1 (in which 
case the top row of the staircase is forced to contain 1, 2, 3). 

The combinatorial interpretation proves in particular that (— l) h A^ > for all 
k > 0. Faulhaber stated this, but he may not have known how to prove it. 

Denoting the determinant by D(w,k), Jacobi's recurrence (**) implies that we have 

(w - kf(w -k + 1)0 - k - l)D(w, k-1) 

= (2w - 2k)(2w -2k-l)(w-k- l)D(w, k) - 2w(2w -l)(w- l)D(w - l,k) ; 

this can also be written in a slightly tidier form, using a special case of the "integer basis" 
polynomials discussed above: 



It does not appear obvious that the determinant satisfies such a recurrence, nor that the 
solution to the recurrence should have integer values when w and k are integers. But, 
identities are not always obvious. 

Generalization to noninteger powers. Recurrence (*) does not require w to be a 
positive integer, and we can in fact solve it in closed form when w = 3/2: 



D[w, k - 1) = TiO - k - l)D{w, k) - T^w - l)D(w -l,k). 
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Therefore A^ 2 ^ = ( 1 { 2 )4 _fc is related to the kth Catalan number. A similar closed form 
exists for ^ m+1 / 2 ) when m is any nonnegative integer. 

For other cases of w, our generating function for involves B n (x) with noninteger 
subscripts. The Bernoulli polynomials can be generalized to a family of functions B z (x), 
for arbitrary z, in several ways; the best generalization for our present purposes seems to 
arise when we define 

B z (x) = x z J2( Z k )*~ kB k, 
k>0 ^ ' 

choosing a suitable branch of the function x z . With this definition we can develop the 
right-hand side of 



J2A^u- k = B 2w ( 

k>0 V 



2 



u~ w 



— — . 2w / , \ — k 

y/1 + AU + 1 \ ^ f 2w\ / y/1 + AU + 1 



ETFF ^ (***) 



as a power series in u 1 as u — > oo. 



The factor outside the s ^S n is rather nice; we have 



y/1 + 4n + 1 \ ^ w (w + j/2\_ j/2 

2 \A / h w+3/2 \ 3 



because the generalized binomial series Bi/ 2 (u 1 ^ 2 ) [5, equation (5.58)] is the solution to 



f(u) 1/2 - f(u)~ 1/2 = u~ 1/2 , /(oo) = l 



namely 



Similarly we find 



/(«) = 



y/1 + 4u + 1 
2^ 





v/1 + 4w + 1 ^ _ ^ -Jfe fj/2-k/2 



3 



3 ~k V J 



u -k/2-j/2 



_^ k_(j/2-k/2-l\_ k/2 _ j/2 
S 2 A J-l 
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So we can indeed expand the right-hand side as a power series with coefficients that are 
polynomials in w. It is actually a power series in -u -1 / 2 , not u; but since the coefficients 
of odd powers of u~ 1 ^ 2 vanish when w is a positive integer, they must be identically zero. 
Sure enough, a check with computer algebra on formal power series yields 1 + A^u- 1 + 
A { 2 w) u~ 2 + A { 3 w) u- 3 + 0(w" 4 ), where the values of A { ™ ] for k < 3 agree perfectly with 
those obtained directly from (*). Therefore this approach allows us to express A^ as a 
polynomial in w, using ordinary Bernoulli number coefficients: 

(( 2w \ B 1 2 V 1 W 1 (k- 1/2- 3-1 

\\2k-l) 2k 1 2 \j) 2k-l-j\2k-l-j-l 

The power series (***) we have used in this successful derivation is actually divergent 
for all u unless 2w is a nonnegative integer, because Bk grows superexponentially while 
the factor 

( 2w \ = r-1 ) k (k-2w-l\ = (-l) k T(k-2w) (~l) fc 2w _! 
\k J { ' \ k J T(k + l)T(-2w) ~ T(-2w) 

does not decrease very rapidly as k — > oo. 

Still, (***) is easily seen to be a valid asymptotic series as u — > oo, because asymptotic 
series multiply like formal power series. This means that, for any positive integer p, we 
have 

fe=0 ^ ' V " / A;=0 

We can now apply these results to obtain sums of noninteger powers, as asymptotic 
series of Faulhaber's type. Suppose, for example, that we are interested in the sum 

n 

#(1/3) = V -— 
k=l 

Euler's summation formula [5, exercise 9.27] tells us that 
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and the parenthesized quantity is what we have called B 2 /s(n + 1). And when u = n 2 + n 
we have B 2/3 (n + 1) = B 2/3 (( a/1 + 4w + l)/2); hence 

^ 1/3) -c(i)~iE4 1/3 V/» 

fc>0 

- 3 1/3 , .5. -2/3 _ _17_ -5/3 , 

— 2 u ^ 36 u 1215 u ^ 

as n — > oo. (We can't claim that this series converges twice as fast as the usual one, 
because both series diverge! But we do get twice as much precision in a fixed number of 
terms.) 

In general, the same argument establishes the asymptotic series 

n 

5> a -c(-«) ~ 4- £4 (a+1)/2 V^)/ 2 - fc , 

fc=i a k>0 

whenever The series on the right is finite when a is a positive odd integer; it is 

convergent (for sufficiently large n) if and only if a is a nonnegative integer. 
The special case a = — 2 has historic interest, so it deserves a special look: 



n 



^f~t-^ 1/2) "" 1/2 -' 4< " /2)m " 3/2 -"' 

k=l 

= - „-!/' + — - i?L U -V a + i2L u -7 /a _ ^1 + . . . . 

6 24 1920 7168 491520 

These coefficients do not seem to have a simple closed form; the prime factorization 32021 = 
11 • 41 • 71 is no doubt just a quirky coincidence. 
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Academia Algebras, a book that is evidently extremely rare: An extensive search of printed 
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at Cambridge University Library, where the volume once owned by Jacobi now resides. 
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at Stanford, so that other interested scholars can take a look.) Ivo Schneider, who is 
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half of Lemma 4, in the case r = 1. And Ira Gessel pointed out that the coefficients in the 
expansion n 2m+1 = J2 a k {^k+i) are cen t ra l factorial numbers in slight disguise. 
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